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Abstract 

The goal of this paper is to identify the universal Whitham hier- 
archy of genus zero with a dispersionless limit of the multi-component 
KP hierarchy. To this end, the multi-component KP hierarchy is 
(re) formulated to depend on several discrete variables called "charges" . 
These discrete variables play the role of lattice coordinates in underly- 
ing Toda field equations. A multi-component version of the so called 
differential Fay identity are derived from the Hirota equations of the 
r-function of this "charged" multi-component KP hierarchy. These 
multi-component differential Fay identities have a well-defined disper- 
sionless limit (the dispersionless Hirota equations). The dispersionless 
Hirota equations turn out to be equivalent to the Hamilton-Jacobi 
equations for the S'-functions of the universal Whitham hierarchy. 
The differential Fay identities themselves are shown to be a generat- 
ing functional expression of auxiliary linear equations for scalar- valued 
wave functions of the multi-component KP hierarchy. 
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1 Introduction 



Nowadays most soliton equations are known to be (re)formulated in the form 
of Hirota equations for r-functions [T] . If such a sohton equation has a dis- 
persionless hmit {dispersionless integrable system), a natural question will 
be whether the Hirota equations have a reasonable dispersionless limit [dis- 
persionless Hirota equations) . This issue has been studied for the following 
relatively small number of cases (because the dispersionless integrable sys- 
tems themselves are, at least currently, rare species). 

1. The first example of dispersionless Hirota equations was discovered for 
the dispersionless KP hierarchy [2j. This dispersionless Hirota equa- 
tion was obtained as a dispersionless limit of the so called differential 
Fay identity [3J (which actually turned out to be equivalent to the KP 
hierarchy itself). Another derivation using a Cauchy kernel was later 
developed [1]. 

2. Dispersionless Hirota equations are also known for the dispersionless 
Toda hierarchy, and applied to Laplacian growth [SI El El E] , WDVV 
associativity equations [9], string field theory [lOl [TH [12l [13] , etc. More- 
over, a connection with some notions of complex analysis (the Grunsky 
coefficients and the Faber polynomials) was pointed out [m [15]. Re- 
cently, a set of Fay-like identities for the Toda hierarchy were derived 
and shown to reduce to the dispersionless Hirota equations in the dis- 
persionless limit [16]. 

3. A few other exotic soliton equation and some problems in complex anal- 
ysis have been studied in the context of dispersionless Hirota equations 
[I7l[T8l[l9l[20]. 

The goal of this paper is to enlarge this list to include a multi-component 
generalization of the KP hierarchy {multi- component KP hierarchy) [211 ESI 
[231 [211 [25]. 

Unlike the aforementioned cases, the multi-component KP hierarchy is 
usually formulated in a matrix form with matrix (pseudo) differential oper- 
ators and vector (or matrix) valued wave functions. This yields technical 
difficulties when one naively attempts to give a prescription of dispersionless 
limit in the Lax formalism. As regards the multi-component KP hierarchy, 
thus, defining a reasonable dispersionless limit itself has been considered a 
tough problem. 

These difficulties can be circumvented in two ways. One way is to resort to 
a scalar Lax formalism based on scalar operators and wave functions [2S] • A 
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prototype of such a scalar Lax formalism can be seen in the Toda hierarchy, 
which is equivalent to the two-component KP hierarchy that depends an 
extra discrete variable s (i.e., the lattice coordinate) [27], and nevertheless 
has a scalar Lax formalism (in terms of difference operators in s). It is rather 
straightforward to generalize this point of view to the case with more than 
two components. The other way is to seek for a prescription of dispersionless 
limit for Hirota equations rather than the Lax formalism. Actually, these 
two approaches lead to the same result. In this paper, we primarily choose 
the route from Hirota equations, which eventually reaches the scalar Lax 
formalism as well. 

In the case of the + 1-component KP hierarchy, we can introduce 
discrete variables Si, . . . ,siy. These variables stand for charges of a state in 
the Fock space of an A^ + 1-component charged free fermion system in the 
field theoretical formalism [231 ES] . In analogy with the aforementioned inter- 
pretation of the Toda hierarchy, this charged multi-component KP hierarchy 
may be thought of as a generalization of the Toda hierarchy as well. 

The dispersionless limit of this hierarchy turns out to be the universal 
Whitham hierarchy of genus zero. The universal Whitham hierarchy was in- 
troduced, at all genera, by Krichever [28] as a universal framework for both 
dispersionless integrable systems and Whitham modulation equations. The 
nonzero-genus cases were recently applied to Laplacian growth of multiply- 
connected domains [221 ED]- The zero-genus case, too, is the target of recent 
studies [211 E21 E3], which have revealed interesting new aspects of this hi- 
erarchy. One can see from our results that many of them stem from the 
multi-component KP hierarchy. 

This paper is organized as follows. Sections 2 is a review of the universal 
Whitham hierarchy of genus zero and some related equations. Of particular 
importance is the Hamilton- Jacobi equations satisfied by the S'-functions. 
These equations play the role of auxiliary linear equations in usual soliton 
equations. In Section 3, we show that these equations can be converted to a 
generating functional form, which we shall later identify with dispersionless 
Hirota equations. Section 4 presents a brief account of the r-function of 
the "charged" multi-component KP hierarchy and the Hirota equations in a 
generating functional bilinear form. In Section 5, we derive a set of differential 
Fay identities for this hierarchy, and show that their dispersionless limit 
(namely, the dispersionless Hirota equations in this case) coincide with the 
equations obtained in Section 3. In Section 6, we show that these differential 
Fay identities are in fact a generating functional expression of auxiliary linear 
equations in the aforementioned scalar Lax formalism. 
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2 Universal Whitham hierarchy of genus zero 



2.1 Dynamical variables 

We consider the universal Whitham hierarchy of genus zero with + 1 
marked points (or "punctures") on a Riemann sphere with coordinate p. 
One of these marked points is fixed at p = oo, and the others are located at 
p = gi, . . . , gAT, which are part of dynamical variables. The other dynamical 
variables are the coefficients of Laurent series Zq{p),Zi{p), . . . ,zj^{p) of the 
form 

oo 

(1) 



oo 

Za{p) = \- Uaj{p - QaY'^ {a = 1, . . . , N) 



defined at the marked points. 

The hierarchy has N sets of time variables ^on = 1, 2, . . .) and i^n 
{a — 1, . . . ,N, n — 0, 1, . . .) attached to the marked points at p = oo and 
p — qa {a. = 1, . . . , N) , respectively. The lowest ones ioi, ^lo, ■ ■ ■ , ^ivo play a 
special role. For convenience, we introduce the auxiliary variable 

N 



too — ~ taOj 



a=l 



which is, of course, not an independent variable. Let dan denote the deriva- 
tives 

9an — 9/Otan- 

2.2 Lax and Zakharov-Shabat equations 

Time evolutions of the hierarchy are generated by the dispersionless Lax 
equations 

danZpip) = {^an{p): Zp{p)} (2) 

with respect to the Poisson bracket 

df dg _ df dg 
^^'^^-^p^to^ dto,dp ^''^ 
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The "Hamiltonians" ^an{p) arc defined as 

^ „^ fn = 1,2,...), 

(4) 



^ao(p) = -l0g(p- ga), 

where ( )(o,>o) denotes the projection to nonnegative powers of p, and 
( )(a,>o) the projection to positive powers of (p — qa)~^- In other words, 

z^{pr = Qan{p) + 0{l) (p^qa) ^ ^ 

for n > 1. In particular, 

^01{P)=P: ^al{p) 



P-Qa 



The Hamiltonians Jlan(p) satisfy the dispersionless Zakharov-Shabat equa- 
tions 

dl3n^am{p) " dam^l3n{p) + {^am{p), ^I3n{p)} = 0. (6) 

They can be converted to the equation 

cuAcu^O (7) 

of the closed 2-form 



oo N oo 

n=l a=l n=0 



By Darboux's theorem, Zf^^p) has a "conjugate" variable C^(p) with which a; 
can be written in the "canonical form" 

u^dzpip) AdCpip). (8) 
Consequently, Cpip) satisfies the Lax equations 

danCpip) = {^an{p), Cpip)} (9) 

of the same form as zpip), alongside the canonical Poisson commutation 
relation 

{zf,ip),Up)} = l- (10) 

Thus (i3{p) amounts to the Orlov-Schulman functions in the dispersionless 
KP hierarchy. 
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2.3 ^-functions 

We can rewrite (E]) as 

d{e + Cpip)dzf,{p)) = o 

where 

oo N oo 



71=1 a=l n=0 

This imphes the existence of a function Sp{p) such that 

ci5^(p) = ^ + C,p{p)dzp{p) 

or, more exphcitly, 

oo N oo 

dS(3 (p) = C/3 (P) C?^/3 (P) + ^On {P) dton + ^^'^^an{p) dtan ■ ( 1 1 ' 

71=1 «=1 n=0 

So{p) and Sj3{p) turn out to have Laurent expansion of the following form: 

°° °° -2o(p)~" 

'5o(p) = X]^On^o(p)" + too\ogZo{p) - ^ -^^ Von, 

71=1 71=1 

n=l 71=1 

where von,Van and 0/3 are functions of the time variables. 
These ^-functions can also be written as 

So{p) = So{zo{p)), Sp{p) = Sf^{zf^{p)) 

where So{z) and S^{z) are defined as 

oo oo _„ 

Z 



(12) 



So{z) = ^ ton + too log ^ - ^ Von, 

n=l n=l 

00 00 _^ 

= ^ t/SnZ"- + tpo log Z + 0/3 - ^ V^n- 



(13) 



n 

n=l 71=1 



The latter functions, too, are called S'-functions, which play a more funda- 
mental role in the subsequent consideration. 



As ffTTj) implies, differentiating ffT2|) by Zj^lp) while leaving the time vari- 
ables constant yields Co(p) and C/3(p)) which thus turns out to be Laurent 
series of the form 



tn 



Up) = Y^ntonzoipr-' + -fr + 5^^o(p)— ' 



n=l 

oo 



n=l 

oo 



(14) 



n=l 



n=l 



Similarly, differentiating (1121) by the time variables gives rise to the following 
expressions of flanipYs: 



Up) - OOnVOm, 



^On{p) 



m=l 

oo 



m 



-m 



a, 



(15) 



m=l 



m 



^an{p) 



-(5„olog2;o(p) - ^ — - — daOVom, 



m=l 



m 



daplogZfS^p) + da0(i)f3 - ^ 



zpipY 



(16) 



m=l 



m 



2.4 Hamilton- Jacobi equations 

Since f2oi(p) = P, the foregoing expression of fioi(p) implies that Zo{p) and 
2;/3(p) satisfy the equations 



(p) 

m=l 

oo 



P = <9oi0/3 - ^ 



m=l 



m 

MpY 

m 



(17) 



In other words, the inverse functions p = po{z) and p = Pf3{z) of 2; = zq{p) 
and z = zp{p) are given explicitly by 



^ — lib 

Po{z) = z-y^ 9oi%m = doiSo{z), 

m=l 

Pfsiz) = doiCpfS - doiVpm = doiSfsiz). 



m=l 
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As a byproduct of these relations, we find an expression of and in terms 
of the coefficients of Sp{z): Since Zfj{p) is assumed to have Laurent expansion 
of the form 

at p = qj3, the inverse function should have Laurent expansion of the form 

pp{z) = qi3 + rpz~^ + 0{z'^). 

Comparing this with ([T]), we readily find that 

<?/3 = ^010/3, = -doivpi. (19) 

We can now derive a system of Hamilton- Jacobi type for the ^-functions. 
This is achieved by substituting p = ppi^z) (/? = 0, 1, . . . , iV) in ffTTl) . The 
outcome of substitution is the equation 

oo N oo 

dSf3{z) = S'p{z)dz + ^ VtQn{Pp{z))dton + X] rian{Pf3{z))dtan, (20) 
n=l a=l n=0 

which implies that 

danSp{z) = naniPfsiz)). 

Since Pfs^z) is related to the S'-functions as (|T8|l shows, we eventually obtain 
the Hamilton- Jacobi equations 

d^nSpiz) = QUdoiSpiz)). (21) 

The lowest Hamilton- Jacobi equations have a special meaning. Whereas 
the equation for tgi is trivial, those for t^o (a = 1, . . . , N) read 

daoSpiz) = - \og{doiSp{z) - g,), (22) 

which can be solved for pp^z) = doiSis^z) as 

Pfsiz) = e-^""^'^^^) + q^. (23) 

In particular, choosing (3 = a and extracting the coefficient of z"'^, we find 
another expression of Tq,: 

r„ = e-^"0'^". (24) 

These Hamilton- Jacobi equations characterize the S'-functions in their 
own terms. One can recover the previous setting of the universal Whitham 
hierarchy by just changing variables from z to p. In this sense, the two 
systems are two different pictures of the same system, and one can move 
from one picture to the other by ([T]) and (ITSll . 
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2.5 F- function 



Following Martinez Alonso, Mafias and Medina [33], we define the F-function 
(the logarithm of the dispersionless r-function [28]) to be a solution of the 
following equations: 

donF = Von, danF = Van (n = 1, 2, . . .), 

aaoi^ = -0a + X^t/3olog(-l) ia = l,...,N), ^^^^ 

(3=1 

where log(— 1) is understood to be equal to, say, ni, though the choice of the 
branch is irrelevant in the final result. Now we can rewrite ffT3D as 



So{z) = J^tonz"" + too\ogz - Do{z)F, 

n=l 
oo 

Sa{z) = ^tanZ"' + t^o log ^ + 0a " Da{z)F, 
n=l 

where Dq{z) and Da{z) denotes the following differential operators: 

Dg{z) = don, Da{z) = (9„„. 



n — ' n 

n=l 71=1 



It should be mentioned that the last part of (l25l) is slightly different from 
that of Mafias et al. [33] ■ Their proof of consistency of the defining equations, 
however, persists to be valid in this form as well. We have modified their 
definition so as to match the fermionic formula of the tau functions by Date 
et al. I23j. 



3 Dispersionless Hirota equations as gener- 
ating functional form of Hamilton-Jacobi 
equations 

3.1 Faber polynomials 

Following Teo's idea |15j developed for the case of the dispersionless KP and 
Toda hierarchies, we now consider the notion of Faber polynomials of po{z) 
and Pa{z), and show that they actually coincide with Qon{p) and Qan{p)- 
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The Faber polynomials $on(p) and $an(p) of Po(-s) and Pa{z) are defined 
by the following generating functions: 



oo 

log ^-^-^ ^ = - > <l>On g , 



n 

n=l 



oo 



(27) 



log = - > $an(?)- 



Rewriting the left hand side as 

.Po{z)-q . po{z) 



log = log h log 1 

z z \ 



Po{z) 
and 

q-Pa{z) Pa{z)-qa 



q-Pa{z 

log '- =log 1 - 

a — \ 



q-qa ^ q-qa 

and recalling that 

Po{z) ^ Z + 0{Z~'^), Pa{z) ^qa + Oiz"'^), 

we readily sec that 

1) ^Qn{q) is a monic polynomial of q, and 

2) ^aniq) is a polynomial of {q — qa)~^ with no constant term. 
Lemma 1. The Faber polynomials coincide with flon{p) and flan{p)- 

$On(g) = ^On{q), $an(g) = ^an{q)- (28) 

Proof. A more explicit expression of these polynomials can be obtained by 
differentiating the generating functions by z as 



Po(z) -q z 



^ oo 



Pa[Z) 



n=l 

oo 



Pa(z) - q ^ 

and extracting the Laurent coefficients by contour integrals. As regards 
^On{z), this yields the contour integral formula 

27n J po{z) - q 2m J\p\=R p-q 
11 



where the contour \p\ = R in the second integral is understood to be suffi- 
ciently large and to encircle q anti-clockwise; the contour in the first integral 
is its image under the mapping p zq{p). The last contour integral is noth- 
ing but the polynomial part of -ZqIp)" evaluated at q, i.e., Qon{q)- In niuch 
the same way, we have the contour integral formula 

^ I z"p'^{z)dz 1 / Za,{pYdp 



2711 J p^[z) - q 2m J\p^q^\=r p-q 

where the contour |p— g^l = r in the second integral is chosen to be sufficiently 
small and to encircle clockwise, leaving q outside. This contour integral 
leaves negative powers of the Laurent expansion of z^ipY at p = ga, hence 
coincides with VLan{,q)- n 
Bearing this interpretation of VtanipY^ in mind, we turn to the Hamilton- 
Jacobi equations ( 12T|) . The goal is to transform these equations, more pre- 
cisely, those for n > 1, into a generating functional form. 

3.2 First subset of Hamilton-Jacobi equations 

Let us first consider the equations for a = 0, /? = and n > 1: 

d,nS^{z) = fion(5oi^o(^)) = ^,n{p,{z)). (29) 

Lemma 2. ( f^) is equivalent to 

^^ rM(z)-rM(w) ^ (30) 

Z — W 

Proof. Substituting q = po{w) in the generating function of the Faber poly- 
nomials ^on{q) = ^onil) yields the identity 

Poiz)-Po{w) ^ z^'" 



log = - > ^ ^On{Po{w)) 



n 

n=l 



Therefore (I29p can be cast into the generating functional form 

PQ{z)-poiw) V^^;-"" 

log = - > donSoiw). 

z ^-^ m 

n=l 

(l26l) implies that dQnSoiw) on the right hand side can be expressed as 

^ — ^ m 
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Thus the generating functional equation can be rewritten as 
Po{z)-po{w) ^ z-"" f „ ^ u;-™ 

n=l m=l 

= log(l~^^ +Do{z)Do{w)F. 

Moving the first term on the right hand side to the left hand side, we obtain 
as a generating functional form of ([211) • ^ 



3.3 Second subset of Hamilton-Jacobi equations 

The second subset of equations are those for a = 1, . . . , N, P = 0, and n > 1: 

danSo{z) = Qan{doiSo{z)) = flaniPo{z)). (31) 

Lemma 3. I[3l\) is equivalent to 

log —-^ = Dq[z)Dp[w)F. 32 

Proof. Substituting q = po{w) in the generating function of the Faber poly- 
nomials ^anil) = ^anil) leads to the identity 

, Po{w)-pa{z) V^^'^o f ( \\ 

log = - > Vt^n[Po{w)). 



log = - > ^ danSo{w). 



Therefore (]3T1) can be converted to the generating functional form 

'o(w) -Pai^ 

Po{w) - g„ 



n=l 

implies that danSo{w) on the right hand side can be rewritten as 



oo 

W 



danSoiw) = - dandomF. 

^ — ^ m 

m=l 

Thus the generating functional equation turns into ( l3Ti) . □ 
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3.4 Third subset of Hamilton-Jacobi equations 

Let us now consider the equations for a = 0, P = 1, . . . , N and n > 1: 

donSpiz) = Qon{doi S 13 (z)) = i^On (P/3 (2) ) • (33) 

Unlike the previous two cases, we now have to use ( !22l) to rewrite these 
Hamihon-Jacobi equations. 

Lemma 4. // f2^) is satisfied, (E^j is equivalent to ^3E) . 

Proof. As in the case of (1291) . we can convert (1331) to the generating functional 
form 



m 

m=l 



log ^ = - > donS,3{w). 

z ^-^ n 

n=l 

By ( l26i) . we can rewrite donSp^w) on the right hand side as 
This leads to an equation of the form 

log P^M^M^ = ± Llg„,a,,F + ± ^^a,,.d,„,F. 

z ^-^ n ^-^ nm 

n=l n,m=l 

The first sum on the right hand side can be evaluated as 

dondpoF = df30 donF 

^-^ n ^-^ n 

n=\ n=l 

= 5/3o(-5'o(2;) +too log 2;) 
= -9/305*0 (^) - logz. 

If ( I22I) is satisfied, we can rewrite dp^^S^iz) as 

dfi^SQ^z) = -log(p/3(z) - qfs) 

and obtain fl32l). □ 
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3.5 Fourth subset of Hamilton-Jacobi equations 

The fourth subset consists of the Hamilton-Jacobi equations for a = (3 = 
1, . . . ,iV and n > 1: 

danSa{z) = ^an{doiSaiz)) = (34) 

Lemma 5. // (2^) are satisfied, [34\ ) is equivalent to 

log 7 P^H-P^(^) ^ ^ Q^^^^ ^ D^{z)D^{w)F (35) 



Proof. As in the case of ( 1311) . we convert (1341) to the generating functional 
form 



log = - > ^ OanSaiw) 

Pa{w) -Qa ~i ^ 

and use ( l26|l to rewrite danSaiw) on the right hand side as 

danSaiw) = - dandaO^ - }] danOamF. 



m 

m=l 



This leads us to an equation of the form 

Pa{w) -Pa{z) 



log 



Pa[W) - qa 

oo 



= logfl - — ) + dandaoF + ^ OandaoF. 

n=l n,m=l 

The first sum on the right hand side can be evaluated as 

oo _^ oo _^ 

dandaoF = OaO ^anF^ 

n=l n=l 

oo 

= <9aO (y-Sa{z) + ^ tanZ"^ + taO log Z + 9^00^) 
n=l 

= -daoSa{z) + log2 + (9„o0a- 

As regards daoSa{z), ( l22l) implies that 

C?aO>S'a(2;) = - \og{pa{z) - qa) ■ 

We now move all logarithmic terms to the left hand side, and obtain ( |35l) as 
a generating functional form of □ 
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3.6 Fifth subset of Hamilton-Jacobi equations 

The fifth subset of equations are those for a, /5 = 1, . . . , A^, a 7^ /5, and n > 1: 

Lemma 6. // [2^) is satisfied, [3^) is equivalent to 

log TT ^ r = 9/3O0a + D^{z)Dp{w)F 37 

We first convert these equations to the generating functional form 

Pp{w)-Pa{z) ^^'""^ \ 
log = - > <9an5/3 W 

Pp{w)-qa ^ « 

and use to rewrite the derivatives on the right hand side as 

^-^ m 

m=l 

The outcome is an equation of the form 

^ ' ra=l n,m=l 

The first sum on the right hand side can be evaluated in the same way as 
the previous case, namely, 

00 _^ 00 _^ 

n=l n=l 

= -df3oSa{z) + <9/3o0O 

= log(pa(2;) - qp) + df30(pa, 

where fl22p has been used in such a form as 

Consequently, we have (1H7|) as a generating functional form of (IHUIl . □ 
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3.7 Dispersionless Hirota equations 

We have thus converted the n 7^ sector of the full Hamilton- Jacobi equa- 
tions ( I2TI) to the four generating functional equations ( 130|) . ( 132|) . ( 135|) and 
( 1371) . Since the n = sector ( l22i) of the Hamilton- Jacobi equations is not 
included therein, these generating functional equations are still incomplete. 

As the following theorem shows, we can modify these generating func- 
tional equations to include as well. The complete generating functional 
equations, ( |38l) - (l4Tl) . resemble the dispersionless Hirota equations for the 
KP and Toda hierarchies. As we shall show later, these equations can be 
derived from bilinear equations of the r-function of the multi-component KP 
hierarchy. For these reasons, we call these equations "dispersionless Hitota 
equations." 



Theorem 1. The full Hamilton- Jacobi equations i \21\) are equivalent to the 
following dispersionless Hirota equations: 

^Do{z)DoMF ^ ^ doi{Do{z) - Do{w))F ^ ^gg^ 

z — w ' 

^gDo(.)(a.o+D.W)F ^ z-doiiDoiz)-d^o~D^{w))F, (39) 
g(a.o+D.(.))(9.o+i5„(«,))F _ zwdoiiD^jz) - D^{w))F ^ 

z — w 

e^^e^9^o+D^iz))id,o+D,M)F ^ -doi{d^o + D^iz)-dpo-D^{w))F,{Al) 
where e^p is a sign factor defined as 

' -1 (« < 13), 



(0) 



-1 {a>/3). 



Proof. Because of the previous lemmas, we have only to show that the gener- 
ating functional form of the Hamilton- Jacobi equations, upon supplemented 
by the n = sector (12^ . is equivalent to the dispersionless Hirota equations. 
As we show below, the four generating functional equations ( !30l) . ( j32l) . ( |35l) 
and ( |371) respectively correspond to the four parts (|38l) . ( l39l) . (HOl) and ( 1411) 
of the dispersionless Hirota equations. 



( I30j) (I38|) To derive flHSl) from flHUl) . we exponentiate the latter as 

^Doiz)Do{w)F 

Po{z) - Po{w) = 1 

z — w 

and substitute 

Po{z) = z - doiDo{z)F, po{w) = z - doiDo{w)F 
(see (|T8l) and (1261) ). This procedure is obviously reversible. 
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(1^ dsn) We exponentiate (fSD as 

and examine the both hand sides. As regards the left hand side, we substitute 

Po{z) = z- doiDo{z)F, pp{w) = -doidf^oF - dQiDp{w)F 
and find that 

LHS = z - doi{DQ{z) - d^o - Dp{w))F. 
As regards the right hand side, we use the relation 

implied by (I23]) and ([26]). Thus we obtain (l39l). 

( I39|) => (I32j) Taking the limit as w — oo in fl39p . we obtain the relation 

z - dm{Do{z) - dpo)F = ze^^"^'^'^^ = e-^^o^°(^). 
Since ( ITQl) and (1251) imply that 

dmdpoF = -qp, po{z) = z- doiDQ{z)F, 
we can rewrite this relation as 

which is equivalent to one of the equations of Having this equation, we 
can recover fl32D from 



(1^ ^ (gO]) We rewrite (ES]) as 

^ = (Pa(^) - <la){Pa{w) - g„) exp((9„o0a + Da{z)Da{w)F) 

z — w 

and examine the both hand sides. The left hand side can be processed in the 
same way as in the preceding cases: 

_ doiiDajz) - D^{w))F 
z — w 

As regards the right hand side, we use the relation 
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and its counterpart for Paiw) — to rewrite it as 

RHS = z-^w-^ exp(-9,o0 + d^oDa{z)F + daoD^{w)F + D^{z)D^{w)F) . 
Since ( 1251) implies that 

dao4>a = -dloF + log(-l), 

we can further rewrite the foregoing expression as 

RHS = -z~'w-^exp{dlQF + d^QD^{z)F + d^oD^{w)F + D^{z)D^{w)F) 
This completes the derivation of ( l40l) . 

(I40h =^ (I35p Taking the limit as w ^ oo in (HOl) leads to the equation 
Since 

doi(pa = Qa, dloF = -daO(p + log(-l), 

we can rewrite the last equation as 

This is equivalent to one of the equations of ( l22l) . We can thereby recover 
(1351). 



( 1371) dUD We first exponentiate as 

Repeating almost the same calculations as the previous cases, we can rewrite 
the left hand side as 

LHS = 9oi(5,o + D^{z) - dpo - Dp{w))F 
As regards the right hand side, we use the relation 

and its counterpart for ppiw) — qa implied by ( l23l) and ( l26l) to find that 
RHS = exp(-9«o0/3 - daod/soF + (d^o + D^iz))id(so + Dpiw))F). 
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On the other hand, fl23|) imphes that 

[0 («>/?), 

so that 

We can thus derive (HT!) from (1371) . 

(14 ip ^ (1371) This is mostly parallel to the previous cases. We take the 
limit as w — > oo in (14T|) . The outcome is the equation 

It is now easy to see that the left hand side is equal to Pa{z) — qp. As regards 
the right hand side, we can use the aforementioned relation 

exp(-9;3o0Q - daodfSQF) = -e„/3 
once again to rewrite it as 

Thus the foregoing equation boils down to 

We can thereby recover flH7|) . □ 

4 Tau function of mutli-component KP hier- 
archy 

4.1 Fermionic formula of tau function 

Following Date et al. [23J (see also the work of Kac and van de Leur [2S]), we 
now consider an A^ + l-component and "charged" version of the KP hierarchy 
in their fermionic formalism. This hierarchy has an iV + 1-tuple t of time 
variables organized in the vector notation as 

^ = (*0; ^1, • • • , ^TV), tcx = ital,ta2, ■ ■ ■) , 
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and extra "charge" variables 

N 

S = {So,Si, . . . ,Sn), y^gg = 0- 

a=0 

The charge variables correspond to the tao's of the universal Whitham hier- 
archy. 

The fermionic formalism is based on an A'^ + 1-component free fermion 

system. The one-particle creation-annihilation operators of this system are 
labeled by a = 0, 1, . . . , as ipajj'ipaj U ^ obey the anticommuation 

relations 

They form an infinite dimensional Clifford algebra. The Fock and dual Fock 
spaces are generated by the vacuum states |0) and (0| that satisfy the anni- 
hilation conditions 

(O|Va,=0 (j>0), V'a,|O)=0 (j<0), 

(o|C, = o (j<o), C,|o) = o (j>o). ^ ^ 

These spaces are decomposed to eigenspaces of the charge operators 

oo 

Hao = ^ : ipaji^aj '■ (normal ordering), 

j=-oo 

each eigenspace being labeled by the charge vector s = (sq, Si, . . . , sjv) men- 
tioned above. We can choose a ground state \s) and its dual (s| in the charge 
s sector as 

{s\ = (0|$o.o*i.i ■ • • ^Ns^, \s) = • • • $L,^o.o|0)' (44) 

where 

'V<o---Cs-i (s>0), 

■ ■'ipa,s (s < 0), 
■ipa,s-l ■ ■ ■ 1pa,0 (s > 0), 

We then have the relation 

{S\lpa,sc,-1 = ^ais){s - Sal, (s|^;_^^_i = (s) (s + 6^ | (45) 
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^aso 



aso 



and a similar relation for \s), where denotes the unit vector 

Bq, = (. . . , 0, 1, 0, . . .) (1 in the a-th component) 
and ea{s) the sign factor 

The r-function r(s, t) of the N + 1-component KP hierarchy is given by 
the expectation value 

T{s,t) = {s\e''^'^g\0) (46) 

of a product of two elements e^^*^ and g of the Clifford groups (i.e., the group 
of invertible elements of the Clifford algebra whose adjoint action generates a 
linear transformation on the linear span of one-particle creation-annihilation 
operators). H{t) is the linear combination 

N oo 



a=0 n=l 

of part of the generators 



Hgn= Yl -M^aj+n- (47) 

jr = -00 

of an + 1-component Heisenberg algebra, g is a general element of the Clif- 
ford group whose adjoint action on the linear span of one-particle operators 
does not mix ipai^ and i'a/s: 

N oo N oo 

a=0 j=—oo a=0 j=—oo 

A typical case is the exponentiated fermion bilinear form 

N oo 



a, (3=0 j,k=—oo 



Note that e^^*\ too, is an operator of this type. 

We have two important consequences of this setting. Firstly, the expec- 
tation value of (H6l) vanishes unless so + si + -- --|-Siv = 0. Secondly, the 
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infinite matrices of tlie coefficients aajjjk and daisjk turn out to be contragra- 
dient to eacli other (i.e., equal to tlie transposed inverse). Tliis implies that 
the operator bilinear identity 

N oo N oo 

7=0 j=— oo 7=0 j=—oo 

holds, and this identity leads to bilinear equations for the r-function. 



4.2 Bilinear equations of tau function 

Let us introduce the free fermion fields 



j=—oo j=—oo 

and rewrite the operator identity fl49l) to an integral form as 
^ f dz ^ f dz 

7=0 7=0 

The contour integral is understood to be an integral along the circle 1^1 = i? 
with sufficiently large radius R: 

- — :Z — On,-l- 

We now apply both hand side of this operator identity to |0) ® |0). A 
clue is the identity 

V7iio)®V';-io) = o, 

which holds for all j G Z because either ip^j or annihilates |0). This 
identity of states in the Fock space implies that 

^ r dz ^ °° 

7=0 7=0 j=— 00 

Thus we are left with the identity 

^ r d7 

Yf^U^)9\^)®r,{z)g\0)=Q. (51) 

7=0 
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We now apply e^^*'-* ® e*^^*^ to the left hand side of the last identity and 
make the inner product with (s' + e^l (S> (s — e^|, where t, t', s, s', a and /3 
are arbitrary. This yields the equation 



E f ^(s' + e„|e^(*'V7(%|0)(^-e;3|e^(*V;(%|0) =0 

7=0 



(52) 



of expectation values. By a multi-component analogue of the well known 
bosonization formula, the expectation values in the integral can be expressed 
in terms of the tau function as 

{s' + e,|e^(*'V7(^)^?|0) = ea^{s')z'''+^"''-'e^^*>^e-''-'^'\{s' + - e^, t'), 
{s + e^|e^(*)^;(z)(7|0) = e^,(s)z-^-+^^--ie-«(*-^)e^-(^V(s - + e„t), 

where D^{z) is the differential operator that has been used in the formula 
(pUj) of the S'-functions, and the other notations are those commonly used for 
the (multi-component) KP hierarchy, namely, 

|-(-l)^.H-i+-+^. (/5>7), ^^^'^ ^ ^" 

Thus the foregoing bilinear equation of expectation values turns into the 
bilinear equation 

^ " dz 



J]e„^(s')e^T,(s) (b 

n ^ 



7=0 ^ 27rz (53) 

X (e'^-("V)(s' + - e^, t'){e^-'^'^T){s - ep + e^, t) = 

of the r-function, which holds for arbitrary value of s', s and for = 
0, 1, . . . , A^. By the standard procedure, this bilinear equation can be con- 
verted to an infinite number of Hirota equations. 

The subsequent consideration is focused on the case where a = /3 = 0, 
which actually covers the general case by shifting s and s' . For convenience, 
we separate the term of 7 = from the sum over 7 = 0, 1, . . . , and move 
the sign factors eoo('S') and eoo(s) to the other terms. €07(5') and €07(5) of 
the latter are thereby replaced by 

~e,{s') = ^ = (-l)<+-+<, ~e,{s) = f^lM = (_i)--^-+^.. 
eoo(s') eoo(s) 
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Bilinear equation fl3^ for a = (3 = thus reads 




(54) 



We shall show later that the sign factor e^(e) is related to the factor propor- 
tional to log(— 1) in the defining equations ( !25l) of the F-function. 

5 Dispersionless Hirota equations as dipser- 
sionless limit of differential Fay identities 

5.1 How to derive differential Fay identity 

Let us recall the differential Fay identity for the r-function r(t), t = (^1,^2, ■ ■ ■), 
of the KP hierarchy [3]. This identity is usually derived by taking a limit of 
the three-term Fay identity 



in the time variables. 

Actually, there is an alternative method to derive the differential Fay 
identity as follows. This alternative method consists of the following steps. 
The first step is to differentiate the bilinear equation 



{zo - Zi){z2 - zs)r{t + [zo] + [zi])r(t + [z^] + [zs]) 
+ (cyclic permutations of 2:1,-22, -23) = 0, 



(55) 



where [z] denotes the shift vector 



[z] = {z,zy2,...,z''/n,... 



/— e«(*'-*'^)(e"''^^V)(i')(e^^^V)(i) = 
J 27r2 



(56) 



by t[. This yields an equation of the form 



;P:e«(*'-*'^)(z(e-^(^)r)(t') + (9ie- 



(^M(0)^(*+r^]) = o, 



where di denotes the differential operator 



di = d/dt 
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The next step is to specialize t' as 

where A and ji are arbitrary constants that sit on the far side of the contour 
1^1 = R oi the integral, namely, 

|A|>i?, \ii\>R. 

The exponential factor e^*^*'~*'^^ thereby reduces to 

g?(t'-t,.) ^ A _ ivVi _ ir' = 

and the foregoing bilinear equation now takes such a form as 

+ (9,r)(* + |A-'| + K'l - \z-']))T(t + \z-']} = 0, 

Assuming, as usual, that the factors including the r-function are holomorphic 
functions of z on the far side of the contour, one can deform the contour 
to a union of circles that encircle the poles at 2; = X, 11,00. The contour 
integrals along those circles are given by residues. Collecting these pieces, 
one eventually obtains the equation 

(A - /.)r(t + [A-i] + ii^-'Mt) - T{t + [X-'Mt + [/x-i]) 

+ (air)(t + [A-i])r(t + [/i-i]) - {d^T){t + [^i-'Mt+[X-']) =0, ^ ' 

which is essentially the differential Fay identity. 

This equation can be cast into a form that is more suited for considering 
the dispersionlcss limit. Let us divide both hand side of this equation by 
(A — /x)r(t + [A"^])r(t + [a*"^]). This leads to an equation of the form 

r(t+ [A-^] + [fi-'Mt] _ ai(logr(t + [A-^]) -logr(f + [/x-^])) 
r(t+[A-i])r(t + [/x-i]) X-^, " ^ ^ 

Using the differential operator 



00 

z 



^(^) = Yl ^n, dn = d/dtn, 



n=l 



we can rewrite the last equation as 



exp((e-W - l)(e-(") - 1) logr(t)) = 1 - Me"'^' - e-^"no,rm 

A — /i 



(59) 



This shows a fully "dispersive" form of the dispersionless Hirota equation 

A — n 
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5.2 Multi-component analogue of differential Fay iden- 
tity 

The foregoing method turns out to work for the bihnear equations fl54p of the 
multi-component well. We thus obtain a multi-component analogue, 

(!6T]) - (!Mll . of the differential Fay identity as follows. 



Theorem 2. The r-function satisfies the multi- component differential Fay 
identities 



(61) 



exp((e^o(^) - l)(e^"('^) - 1) logr(s,t)) 

^ doi (e^°(^) - e^«(^)) log r(s, t 
\ — fi 

Aexp((e^o(^) - l){e^'^o+D^M _ i)iogr(s,t)) 

= A - aoi(e^°(^) - e^"°+^"(^))logr(s,t), (62) 
exp((e^-o+i).(A) _ i)(ga.o+o4M) _ i)logr(s,t)) 

X^idoi [e^'^o+D^W _ ^d^o+D^M^^ log r(s, t) 

= T ) (63) 

A — jj, 

e„;3exp((e^-°+^'^(^) - l){e^^o+D,M _ iogr(s,t)) 

= -9oi(e^"°+^°(^) - e^'3«+^''('^)) logr(s,t), (64) 

where e^"° is understood to be the shift operator in the s variables, i.e., 

e^-^f{s)= f{s-e, + e^). 
Proof. To derive flUTj) . we differentiate (13^ by t'^^ and put 

t'o = to+M + K'], t'^ = tc. (a = l,...,iV), s' = s. 
The differentiated bilinear equation thereby simplifies as 
dz Xfi 



+ (ft,iT)(s, t + |A-']„ + |^-i]„ - |2-']o))r(s, t + lz-%) = 0, 

where [2;]^ now stands for the shift vector 

[z]a = (. . . , 0, [2;], 0, . . .) {[z] in the a-th component) 

in the multi-component time variables. This is substantially the same bilinear 
equation as mentioned above for the case of the usual KP hierarchy. We thus 
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obtain the equation 



r(g, t+[A-^]o + [/x-^]o)r(g,t) 

_ doi{\ogT{s,t+[\-%)-\ogT{s,t+[fi-%)) 



which we can readily rewrite to ( 16T1) . 

To derive ( l62l) . we differentiate ( 15^ by and put 



The differentiated bihnear equation talces such a form as 



+ (9oir)(s - Co + e„, t + [A - [z % + x 
X r(s, t + [2;-i]o) 



X r(s - Co + Bq,, t + [2; = 0. 
After calculating the contour integrals, this equation reduces to 

Ar(s - eo + Ba, t+ + [Ai"^]a)r(s, t) 

- Ar(s -eo + ea,t+ [i2-\)t{s, t + [A"^]o) 

- {doiT){s - Co + e„, t + [/i"^]^)r(s, t + [A"^]o) 

+ (9oir)(s, t + [A"^]o)r(s - eo + e„, t + [/x-^]^) = 0, 

which we can further rewrite as 

t{s - Co + e„, t + [A-^]o + [i^~%)t{s, t) 



= A - 9oi(logr(s, t + [A"^]o) - log r(s - Bq + e^, t + [fi'^]a))- 



After some more algebra, this equation becomes (lU^ . 
To derive ( l63l) . we choose 




(2r(s - eo + e^, t + [A ^Jo - [2; + [/U + 




So -2, s^ = Sc, + 2, 4 = ■^/^ (/^7^0'«)- 
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The differentiated bilinear equation now reads 

/dz 
— (2;r(s - 2eo + 2e,, t - [z-% + [A'^], + + 

+ (9oir)(s - Co + e„, t - [z-\ + [\-\ + x 
X T{s,t+[z~%) 

X r(s - eo + e„, t + [z~^]o) = 0. 
The outcome of residue calculus is the equation 

t{s - 2eo + 2e^, t + [X-% + [fi-%)T{s, t) ^ A/i ^ 
r(s - Co + e^, t + [A-i]a)r(s - Bq + e^, t + A - /i 

X 9oi (log r(s - Co + e^, t + [A"^]q) - log r{s - eo + e^, t + , 

which becomes (1631) after some algebra. 

Lastly, we derive flMl) . This is achieved by choosing 

*L = *« + M, t;5 = */3 + *; = *7 (7 7^«,/3), 

s'o = So - 2, = + 1, s'p = Sp + 1, = (77^ tt, 

and repeating the same procedure as the previous cases. Unlike those cases, 
however, we now have to account for the sign factors in (1541) . The outcome 
takes a different form for a < (3 and a > (3 as follows: 

dz 1 



27ri z^ ^^'^^^ " + f^a + ep,t- [z % + [A % + [fx + 

+ {doiT){s - 2eo + e„ + e^, t - [z^^ + [X'% + [fi~%)) x 
X T{s,t+[z-\) 

/dz 1 A 
77- r(^oi7-)(s - eo + e^, i + [A"^]^ - [z'^]a + [/U^^]/?) x 
ZTTi 2 ^ — A 

X r(s - Co + Cc,, t + [2;~^]a) 

77^ ^(9oir)(e + e^,, t + [A^^]a + [/i"^]/3 - x 

zm z z — fi 

X r(s - Co + e^, t + [z-^]f3) = 0, 

where 6^/3 stands for the same sign factor as in (jHj). This equation reduces 
to 

^ t{s - 2eo + eg + e^, t + [A~% + [/z~^]^)r(5, t) 
°^r(s - Co + e^, t + [A-i]c,)r(s - + e^, t + [fi-^]^) 
= -doi (log r(s - Co + e„, t + [A"^]«) - log r(s - cq + e^, t + , 

and we eventually obtain (IMl) . □ 
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5.3 Dispersionless limit 

Dispersionless limit is realized as quasi-classical limit. Namely, we allow the 
r-function to depend on the "Planck constant" h as well, i.e., r = T{h, s, t), 
and assume that the rescaled r-function 

Tfi(s, t) = T{h, fr^s, fr^t) 

behaves as 

Tn{s,t) = exp{h-^F{s,t) +0{h-^)) (65) 

in the classical limit h 0. 

The rescaled r-function satisfies (j6Ti) - (!64l) with the operators e'^"^^-' etc. 
being also recalled: 

exp((e^^°(^) - l)(e^^°(^) - 1) \ogTn{s,t)) 

^ hdoi (e'^^oix) _ e?iA)(M)) log Tn{s, t) 
X — fi 

= A - ;i9oi(e^°(^) - e'^-o+'^^"('^)) logr(s,t), 

A/in9oi(e^"»+^^"(^) - e^-o+^^"(^)) logra(s,t) 
A — /i ' 
e„^exp((e^"«+^^"(^) - l){e'^i3o+t^D^M - 1) logr;i(s,t)) 

= -^^01 (e'^'^°+^^'^(^) - e^^''+''^^^'''>) logTn{s,t). 

If we substitute the aforementioned quasi-classical ansatz of the rescaled r- 
function and take the limit as /?. ^ 0, we end up with the dipsersionless Hirota 
equations fl38l) - fHTl) for the F-function. In other words, the multi-component 
analogues flUT|) - flM|) of the differential Fay identity are "dispersive" counter- 
parts of the dispersionless Hirota equations. 

6 Differential Fay identities as generating func- 
tional form of auxiliary linear equations 

6.1 How to derive auxiliary linear equations 

It is known that the differential Fay identity is actually equivalent to the 
KP hierarchy itself. This is a consequence of the fact that the auxiliary 
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linear equations of the KP hierarchy can be derived from the differential Fay 
identity. Let us briefly review this fact 0, Appendix C] (see also Teo's paper 
[16] for an analogous result and its proof for the Toda hierarchy). 

To derive auxiliary linear equations, we write the differential Fay identity 

as 

T{t + [A-i])r(t + [/.-I]) A - /.''^ T{t + [A-i]) 

and shift the variables as i ^ i — [A^^] — [/x^^]. This leads to another form 

r(t-[A-V[/x-^])r(t) 1 , r(t-[A-^]) 

r(t - [A-i])r(t - [/i-i]) A - /.''^ r(t - [/.-i]) 

of the differential Fay identity. We can rewrite it as 

^,,, r{t-\^,-^]) r(t-[/.-^]) / 1 ^ r(t-[A-^]) 



r(t) r{t) X-r ^^r{t-[^,-^]) 

which turns into the equation 

Ae-^(^)*(t,/^) = -9ilog^(t,A))^(t,/i) (66) 
for the wave function 

nt,z) = "-^^^e^'^-\ (67) 

fl66|) may be thought of as a generating functional form of auxiliary linear 
equations for ^!{t,z). Let us introduce the fundamental Schur functions 
hj{t), J = 0, 1, 2, . . ., by the generating function 

oo 
j=Q 

(jUUj) can be thereby decomposed to the linear equations 

h,{-dt)^{t,fi)=f^{t)^{t,f,) (j=2,3,...), (68) 
where dt denotes the vector 

\ 2 n 



31 



of derivative operators, and /j (t)'s are the coefficients of Laurent expansion 
of di log z) ai z = oo, i.e., 

oo 

d, logvl/(t, z) = z + 9i(e-^(^) - 1) logr(t) = z + J2 

i=i 

It is well known that these Laurent coefficients give conserved densities of 
the KP hierarchy. 

We can convert fl66l) to linear equations of the familiar evolutionary form 
dn'fit,z) = Bn^it,z) (n = l,2,...), (69) 
where Bn is a differential operator of the form 

Bn = + 6„25r' + ■ ■ ■ + bnn. 

This is achieved by the following method [2, Appendix C]. Let us rewrite 
(|66|) as 

(l-e-^W)vI/(t,/i)=X(t,A)vl>(t,/i), 

where X(A) denotes the differential operator 

X{t, A) = X~^di + 1 - X~^di log ^(t. A) 

that depends on A. Let Xk(t, A), = 2, 3, . . ., be differential operators defined 
by the recurrence relation 

Xfc+i(t, A) = Xk{t, A) - Xk{t - [A-i], A)(l - Xk{t, A)), Xi(t, A) = X{t, A). 
Xfc(t, A) is a differential operator of the form 

Xfe(t,A) = A-'=5f + ■■■ , 
and turns out to satisfy the equation 

(l_e-^5(A))fcvl/(t,^)=Xfc(t,A)vl/(t,/x). 
Moreover, by construction, 

Xfc(A) = 0{X-'). 

Consequently, the action of 

DiX) = - log(l - (1 - e-^W)) = E - ^^''^'^)' 

ifc=i 

on ^(t, /i) is well defined, and we have the linear equation 

oo ^ 

D(A)vl/(i,/i) = E-Xfc(t,A)vI/(t,/i), 

k=l 

which, upon expanded in powers of A, gives the usual auxiliary linear equa- 
tions (l69l) as expected. 



32 



6.2 Generating functional form of auxiliary linear equa- 
tion for mult i- component case 

We now turn to the the muhi-component case. Let us recall the following 
form of the differential Fay identities first derived from the bilinear equations 



r(s, t+[A-i]o)r(s,t+[/i-i; 







1 + T ^oilog^ — . , r. n V 70 

A-/i r(s, t+ A-i o) 



t{s - ep + eg, t + [A + [/i ^]a)r(s, t) 
r(s, t + [A-i]o)r(s - Bq + e^,, t + 



A + 5oilog . 4 I r\-ii ^ ' 

r(s, t + A ^ o) 



r(g - 2eo + 2e^, f + [A'^ + [i2-\)t{s, t) 
r(s - eo + e^, t + [\-'^]c,)t{s - + e^, t + [/i-i] 



A - r(s - Co + Bq,, e + [A'^ 

t{s - 2eo + + e^, t + [\'\ + [/i~^]/j)r(e, t) 
r(s - eo + Bo, t + [A-i]a)r(s - Cq + e^j, t + [/i-^]/?) 



= ao.log -i'--° + -^-*+if:'^j . (73) 

r(s - eo + e„, t + [A i]^,) 

We now convert these equations to auxiliary linear equations for the scalar- 
valued wave functions 

(74) 

To simplify notations, these wave functions are referred to as '^q{z) and 
^^(z). We have the following analogue of (lUBl) for these wave functions. 

Theorem 3. ([7^-([75|j can he converted to the linear equations 

Ae-^(^)^(/i) = -(9oi-9oilog^o(A))^(/i), (75) 
g-a.o-D.(A)^(^) = (9oi-9oilog$,(A))^(/x) (« = l,...,iV) (76) 

for ^(/i) = ^o(/i), ^i(/^), • • • , ^^N{^^)■ 
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Proof. Let us first consider flTUl) . This equation has the same structure as 
the differential Fay identity of the KP hierarchy. Therefore we can repeat 
the foregoing calculations for the KP hierarchy to derive the equation 

Ae-^°W^o(/i) = -{doi - 9oi^o(A))^o(/i), 

which is nothing but (17^ for = "^oin)- 

As regards ( ffTl) . we shift the variables as 

S ^ S + So - Ba, t ^ t - [\^^] - [^~^] 

and consider the equation 

r(s + Co - Cq, t - [A"^]o - [/i~%)r(s, t) 



A 



r(s, t - [A-i]o)r(s + Cq - e^, t - [/i-i]^) 
T{s,t-[X-%) 



A + 9oi log ■ 



thus obtained. We can rewrite this equation in two different ways as 

T{s,t) 

r(s, t - [A"^]„) / r(s + Co - So, t - [A~^ 



and 



A + (9oi log — - — - — 



0, 



t{s + Bo - Ba, t - f r(s + eo - e„, t - [A~^]o' 



(|a + doi log 



r(s,t) V r(s, t - 

Upon interchanging A and /i, the first equation reduces to 
g-9.o-i?.(A)^^^^) = (9oi -9oilogM/„(A))M/o(/i), 

which coincides with fl76l) for the case of ^^(/u) = \E'o(/^)- On the other hand, 
the second equation becomes 

Ae-^«W^„(/x) = -{doi -doi log^o(A))^a(^), 

yielding fl73|) for the case of = 

We now turn to (1721) . By shifting the variables as 



s + 2eo - 2eQ, t -> t - [A % - [/i 
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we obtain the equation 

t{s + 2eo - 2e^, t - [A'^ - [/x-%)r(s, t) 

t{s + eo - e„, t - [X-%)t{s + bq - e^, t - [/i-^]a) 

A/i r(s + Co - t - 

- -doi log 



A-/X r(s + eo - e«, t - [/x-i]a)' 

equivalent ly, 

ris,t) 

_ t{s + Bq - Ba, t - [fi~%) A/i r(s + eo - e^, t - [A"^]a) 

-Coi log ■ 



r(s,t) A-/i r(s + Co - Cq, t - [// 

We now use the difference relation 

e«(s + eo - e^) = -e„(s) 
to convert the last equation to the linear equation 

g-a.o~D.(A)^^(^) = (^01 -9oilogvI/,(A))vl/,(/i) 

for \&Q,(/i). This is exactly ( 1761) for the case of ^^(/i) = "^aifJ^)- 
Lastly, we consider fl73|) . Upon shifting variables as 

s ^ s + 2eo - - 6/3, t^t - - [/i"^]/3, 

we have the equation 

t{s + 2Bo -B^-Bis,t- [A"l]« - [lJ-^]f3)T{s, t) 



t{s + Bo - Be,, t - [X-%)t{s + Bo - B/s, t - [/i"^]/?) 

r(s + Bo - Ba, t - [X'%) 



= <9oi log ■ , . . , 

equivalent ly, 

-a.o-D„(A) + 60 - 6/;, t - [/i-^];3) 



T S,t) 



^ r(g + 6o - 6/3, t - [/i ^]f3) ^ t{s + Bo-Bc,,t- [A %) 

r(s,t) °^ r(s + 6o - 6/3, t - [/i-i]/3) ■ 

Noting the difference relation 

6/3(5 + 60 - Ba) = -eafsifsis), 

we can convert the last equation to the linear equation 

e-^"«-^"(")^/3(/i) = (5oi -9oi log^,(A))^/3(/i) 

for ^/3(/i). This gives (!7U|) for \I'(/i) = ^E'/3(/i), /5 7^ a. Thus all equations of 
(1751) and ( |76l) have been derived from the multi-component differential Fay 
identities. □ 
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6.3 Evolutionary form of auxiliary linear equations 



We can convert the generating functional auxiliary linear equations ( !75l) - 
(1761) to an "evolutionary" form. As it turns out below, those auxiliary linear 
equations turn out to be a mixture of those of one-component KP hierarchies 
and two-dimensional Toda field equations. The KP hierarchies live in each 
sector of the -|- 1 sets of time variables tQ,ti, . . . ,ti\f. The Toda field 
equations (and, actually, the Toda hierarchies) connect the tg-sector with 
the other N sectors pairwise. We can thus eventually recover all building 
blocks of the scalar Lax formalism of the multi-component KP hierarchy 

ESI. 



6.3.1 KP hierarchy in to-sector 

Let us note that the first set ( 1751) of these equations have the same structure as 
their counterpart (IMl) in the KP hierarchy. Therefore they can be converted 
to linear equations of the form 

don^iz) = Bon{doi)^{z) (77) 

for "^{z) = "^oi^z), \E'i(2;), . . . , \E'7v(-2)- -Son(<9oi) is a differential operator of the 
form 



Bon{doi) = di;, + O{dl 



■or') 



Among the A^ -|- 1 wave functions, "^oiz) plays the role of the wave function 
in the usual KP hierarchy. If we define the dressing operator 



01 



such that 



oo 

^o(^) = Wo{doi)z'°e^^'°''^ = (l + Y,^o, 



the auxiliary linear equations for \E'o(^) can be converted to the so called Sato 
equations 

^^^^ = Bon{doi)Wo{doi) - Wo{doi)d^,. (78) 

As usual, this implies that i?on(c^oi) is given by the KP-like formula 

5o„(9oi) = (Lo(9oi)")>o, (79) 
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where ( )>o denotes the projection to nonnegative powers of doi, and Lo{doi) 
is the pseudo-differential Lax operator 

Lo{doi) = Wo{doi)-doi-Wo{doi)-\ 
6.3.2 Toda field equation in {to, tQ,)-sector 

We now turn to ( 1761) . As we show below, the auxiliary linear equations of 
Toda fields emerge in the lowest and next-to-lowest orders of expansion in 
powers of A. 

The Toda fields (pjs = 4>i3{s, t) are defined in an exponentiated form as 
which is actually the leading coefficients of the amplitude part of p{^z\ i.e., 

oo 

i=i 

Let us examine the lowest orders of expansion of (1761) in powers of A. The 
A^'-terms yield a linear equation of the form 

(/X has been replaced by 2;.) Introducing a new field as 

ga = ^Ol'/'a, (81) 

we can rewrite this equation as 

aoi^;3(^) = (e-^"«+ga)*/3(^). (82) 
Let us note that the associated Hamilton- Jacobi equation 



in quasi-classical approximation is nothing but ( l22l) . The A^-terms of ( 1761) 
give the somewhat ugly equation 

-e-^°«9«i^^(2;) = doidai logr(s + Cq - e„, t) ■ ^^{z). 

We can rewrite it as 

d^i^^iz) = ry-^'^H^iz), (83) 
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where 

ra = -doidailogT{s,t). (84) 
The associated Hamilton- Jacobi equation reads 

Using the foregoing Hamilton- Jacobi equation for the toi flow, we can rewrite 
it as 

T 



This coincides with fl?Il) for n = 1. 

As one can thus see from quasi-classical approximation, the auxiliary 
fields Qa and correspond to those of the universal Whitham hierarchy 
defined by ( |T9l) . Moreover, substituting 

in f lHHj) for the case oi (3 = a and extracting the coefficient of ^r^^e^*^*"'^-*, we 
obtain another expression 

r„ = exp(0c.(s, t) - (f)ais - eo + e«, t)) (85) 



of r^. This expression corresponds to ([2] 

f l82|) and fl83l) may be thought of as auxiliary hnear equations of the 
two-dimensional Toda fields with continuous variables toi^^ai and discrete 
variable Sa- The zero-curvature equation 

[doi - e-^"° - g„ d^i - r«e^"°] = 

comprises the two equations 

9oir„(s, t) - {qa{s, t) -qa{s-eo + e^, t))ra{s, t) = 0, ^^^^ 
daiQais, t) + r„(s + eo- e^, t) - rc,(s, t) = 0. 

The first equation of 0861) is automatically satisfied under the definition of 
Qa and Tq, by (IHTl) and (IH31) : the second equation reduces to the Toda field 
equation 

doidai(pais,t) +exp(0c,(s + eo - e^, t) - 0a(s,t)) ^^^^ 
- exp(0„(s,t) - (f)a{s - Co + Sa, t)) = 0. 

This explains why a multi-dimensional dispersionless Toda field equation (or 
the Boyer-Finley equation) shows up in the universal Whitham hierarchy 
[32]. 
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6.3.3 Auxiliary linear equation of Schrodinger type in (tg, ta)-sector 

The (to, t«)-sector is also accompanied by an auxiliary linear equation of the 
form 

{{doi - qa)dai - ra)^(^) = (88) 

satisfied by ^{z) = '^o{z),'^i{z), . . . ,'^]\i{z). These are analogues of two- 
dimensional "integrable Schrodinger equations" [Ml 135]. 

( IHHj) can be derived from the Toda-like auxiliary linear equations 0821] and 
(!H^ as follows. We first apply the operator doi — qa to both hand sides of 
(!83l) . This yields an equation of the form 

{doi - qa)'^p{z) = (9oi - qa){ry-'^p{z)) 

= (9oir,)e^"«^;3(^) + r,(9oi - g,)(e^"»^;3(^)). 

By the first equation of 0861) . we can rewrite the first term in the last line as 

(9oir„)e^"«^^(2;) = (g„(s, t) - g„(s - Cq + e„, t))r^e^-''^ p{z). 

As regards the second term, we use 0821) to rewrite its main part as 

{doi - g„)(e^'^''^^(z)) = e^-'do,^p{z) - g,e^-^^(z) 

= e^'^o(e-^'^o + g„)^;3-ga^/3(^) 

= -^f^iz) + {q^{s - eo + e„, t) - q^{s, t))e^"°*^(z). 

Thus the terms proportional to ga(s, t) — ^^(s — Co + eQ,, t) cancel each other, 
and (IHH]) follows. 

6.3.4 KP hierarchies in other sectors 

Let us now consider 0761) . Applying the shift operator e^"" to both hand side 
yields 

e-^.{A)^(^) = (5oi-e^"°9oi log^„(A))e^"»^(/i). 

On the other hand, (!83l) implies that 

e^"°^(/i) =r-ia„i^(/i). 

Therefore we can eliminate the difference term e^"''\E'(/i) and obtain the dif- 
ferential equation 

e-^4A)^(^) = (^01 - e^-9oi log^,(A))r-i9,i^(/i). (89) 
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Letting as A ^ oo, we have 

^(/i) = (9oi-e^-°g,)r-Ui^(/^). 

Taking the difference of these two equations and doing some algebra, we 
eventually obtain an equation of the form 

(l_e-^'^W)vl/(;,)=X„(A)vI/(^), (90) 

where Xa{X) is the first-order differential operator 

A = 9oilog ■ r„ dai. 

V r(s - eo + e„, t) / " 



We can now repeat the reasoning in the derivation of ( 1691) from (1661) . The 
outcome is a system of "evolutionary" linear equations of the form 

d^^^iz) = B^^{d^,)^iz), (91) 

where Ban{dai) is an n-th order differential operator without 0-th order term, 
i.e., 

ann—l 

The absence of 0-th order term is a consequence of the structure of Xa{z). 
( !69l) may be interpreted as auxiliary linear equations of the one-component 
KP hierarchy in the t^-sector. 



7 Conclusion 

We have thus identified the universal Whitham hierarchy of genus zero with 
the dispersionless limit of the charged multi-component KP hierarchy. In 
the course of this rather lengthy consideration, we have seen a number of 
remarkable aspects of these systems in themselves. 

A main conclusion of this consideration is that the dispersionless Hirota 
equations and the differential Fay identities can play the role of master equa- 
tions that characterize the (usual and dispersionless) r-functions of these in- 
tegrable systems. This interpretation is complementary to the conventional 
point of view based on the Hamilton- Jacobi equations and the auxiliary linear 
equations. We have an honest impression that the approach from the Hirota 
equations is technically simpler and conceptually more essential. This ap- 
proach deserves to be pursued further. 
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Alongside these fundamental aspects, we expect some applications of the 
results of this paper. For instance, the recent paper [13] of Bonora et al. 
conjectures the relevance of a multi-component version of the disperesionless 
Toda hierarchy in light-cone string field theory. This conjecture might be 
explained in our framework. Another intriguing issue is to generalize the 
notion of associativity equations ^ to the mult i- component setting. In this 
respect, we should mention the recent work of Konopelchenko and Magri [SB] , 
in which a novel approach to the universal Whitham hierarchy is discussed 
along with some other issues on dispersionless integrable systems. 
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